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Abstract. From the starting point of the wcU known Reynolds number of fluid 
turbulence we propose a control parameter R for a wider class of systems including 
avalanche models that show Self Organized Criticality (SOC) and ecosystems. R is 
related to the driving and dissipation rates and from similarity analysis we obtain a 
relationship R ^ N^" where N is the number of degrees of freedom. The value of 
the exponent 13m is determined by detailed phenomenology but its sign follows from 
our similarity analysis. For SOC, R = h/e and we show that /5jv < hence we show 
independent of the details that the transition to SOC is when i? ^ 0, in contrast 
to fluid turbulence, formalizing the relationship between turbulence (since [3^ > 0, 
R — > oo) and SOC {R = h/e 0). A corollary is that SOC phenomenology, that is, 
power law scaling of avalanches, can persist for finite R with unchanged exponent if 
the system supports a sufficiently large range of lengthscales; necessary for SOC to be 
a candidate for physical systems. We propose a conceptual model ecosystem where 
R is an observable parameter which depends on the rate of throughput of biomass or 
energy; we show this has Pn > 0, so that increasing R increases the abundance of 
species, pointing to a critical value for species 'explosion'. 
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1. Introduction 

A central idea in physics is that complex, often intractable, behavior can be quantified 
by a few measurable control parameters. In fiuid turbulence a single control parameter, 
the Reynolds number Re quantifies the transition from ordered (laminar) to disordered 
(turbulent) flow. This parameter is from dimensional arguments a function of 
macroscopic system variables, but is also expressible as a function of the number 
of energy carrying modes or degrees of freedom (d.o.f.)[T]. Although the detailed 
dynamics are different, many of the macroscopic features of idealized, finite (large) 
Reynolds number turbulence are also characteristic of other strongly correlated, out 
of equilibrium systems. These systems can all be driven into a disordered state with 
defining characteristics: they have many degrees of freedom (d.o.f.); are driven and 
dissipating, are out of equilibrium but on average in a steady state, and show anomalous 
scaling over a large dynamic range. Loosely speaking, a 'class' of such systems will show 
an order- disorder transition, captured by varying a single control parameter. This 
includes avalanche models exhibiting Self Organized Criticality (S0C)[2l [3l HI O [6l [7] . 
It was originally argued[2l El] that these systems self organize to an SOC state which in 
the above sense is their state of maximal disorder. Subsequent analysis has established 
a consensus 1 1 H IT2 | [T3l |6| [5] that SOC is a limiting behavior in the driving rate h and the 
dissipation rate e, such that h/e with h,e 0, (and h < e, that is, steady state). 
This is exemplified by the constructive definition (in[4j) as "slowly driven interaction 
dominated thresholded" (SDIDT) systems. 

In this paper we give a prescription for obtaining a control parameter for these 
systems, in analogy to the Reynolds number in fiuid turbulence. For avalanche models 
exhibiting SOC, we identify d.o.f. with realizable avalanche sizes and we show that the 
relevant control parameter is h/e. It follows that the SDIDT limit is reached by 
taking Ra to zero; we show this maximizes the number of d.o.f. in the opposite sense to 
fiuid turbulence. This result clarifies the much debated relationship between turbulence 
and SOC[14l[l5l|16l[l7l|lHl[l9]. An important corollary is that SOC phenomenology can 
quite generally persist under conditions of finite drive in a sufficiently large bandwidth 
system. As our result flows from dimensional analysis it is quite generic. An important 
example that we give here IS ClS db parametrization of ecosystem models for species 
abundance [HI O |10]. For ecosystems we show that as the control parameter increases so 
does the abundance of species, or d.o.f. This points to the possibility of a critical value 
at which the onset of diversification of species occurs. 

2. Similarity analysis and Reynolds number 

The systems that we have in mind all have strongly coupled d.o.f. that transport 
some quantity from the driving to the dissipation scale. Provided that this 'dynamical 
quantity' is governed by a conservation law to ensure steady state (not necessarily 
equilibrium) on the average we insist that its precise nature, and the microscopic 
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Table 1. 11 theorem applied to homogeneous turbulence. 



Variable 


dimension 


description 


Lo 


L 


driving length scale 


V 


L 


dissipation length scale 


U 


LT-^ 


bulk (driving) flow speed 


V 




viscosity 



details of how it is transported are not relevant to the macroscopic ensemble average 
behavior. We seek a control parameter expressible in terms of the number of d.o.f. of 
the system that parameterizes the transition from ordered (few d.o.f.) to disordered 
(many d.o.f.) behavior. We identify the control parameters of the system in terms 
of known macroscopic variables by formal dimensional analysis (similarity analysis 
or Buckingham 11 theorem, see e.g. [20]). Any system's behavior is captured by a 
general function F which only depends on the relevant variables Qi„v that describe the 
system. Since F is dimensionless it must be a function of the possible dimensionless 
groupings I^i..m{Qi..v) which can be formed from the Qi.y. The (unknown) function 
F(ni, 112, --^m) is universal, describing all systems that depend on the Qi„v through the 
'^i..m{Qi..v) and the relationships between them. If one then has additional information 
about the system, such as a conserved quantity, the ^i„m{.Qi..v) can be related to each 
other to make F explicit. Thus this method can lead to information about the solution of 
a class of systems where the governing equations are unavailable or intractable, often the 
case for complex systems where there are a large number [N here) of strongly coupled 
d.o.f.. If the V variables are expressed in W dimensions (i.e. mass, length, time) then 
there are M = V — W dimensionless groupings. 

Here, since we have that the precise nature of the transported dynamical quantity 
is irrelevant, the only relevant dimensions are length and time so that W = 2. We next 
insist that there is a single control parameter {R, in the case of turbulence, the Reynolds 
number Re) which may be expressed as a function of the number of active degrees of 
freedom A^. This means that the system's behaviour is captured by some F(ni,n2); 
where R = Hi and 112 = f{N) and the Hi and 112 are related via some conservation 
property. Hence M = 2 so that V = 4; there are four relevant variables to consider. 

It is useful to fix ideas in terms of a relatively well understood example of the above, 
namely turbulence. Our aim here is to obtain a control parameter R by analogy to Re 
via dimensional analysis; for a detailed discussion of the universal scaling properties of 
Kolmogorov (K-41) turbulence and their origin in the Navier Stokes equations see for 
example [1]. As above, for K-41 we have four relevant macroscopic variables (given in 
Table 1) and two dimensionless groups: 

Ui = ^ = RE,U, = ^ = fiN) (1) 

Hi is the Reynolds number of the flow, and the ratio of lengthscales H2 is directly 
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related to the number of d.o.f. available and we now relate Re to f{N) (or Hi to 
112). For incompressible fluid turbulence, our dynamical quantity is the time rate of 
energy transfer per unit mass Si through length scale /. Conservation and steady state 
imply that in an ensemble averaged sense this is balanced by the rate at which energy 
is transferred to the fluid Emj ~ U^/Lq which in turn is balanced by the dissipation of 
energy within the fluid Ediss so that Sinj ~ ~ £diss- Dimensional arguments (e.g. [1]) 
lead to Ediss ~ i^^/rj^ ■ Conservation, that is Einj ~ Ediss then gives the well known result 
[1] which relates Hi and 112: 

The 4/3 exponent that arises for K-41 is modifled if we consider other turbulent 
flows with different phenomenologies, for example anisotropy, and intermittency. 
Nevertheless, for any turbulent flow we can anticipate that the relationship between 
Lo/t] and the number of degrees of freedom N will be of the form: 

N.[^J (3) 

with a > 0; the crucial point is that for turbulence, N always grows with Lo/rj. The 
only property of turbulence with which we are concerned here is that 

Rj^ ^ (^^j (4) 

and that in particular, for turbulence = Pet > 0. This identifles the Reynolds 
number as the control parameter for a process (turbulence) which simply grows more 
active modes or d.o.f. as we increase Re, taking the system from order (few d.o.f. or 
laminar flow) to disorder (many coupled d.o.f.). 

We will now see that more generally, similarity analysis is sufficient to obtain the 
relationship between the control parameter R and the number of degrees of freedom N 
of the form: 

R ~ AT^^ (5) 

The value of the exponent jS^ will depend on the details of these systems but crucially 
we will see that the sign of /^at is obtained from similarity analysis. This is sufficient to 
establish if, as in the case of turbulence, increasing R increases the disorder or complexity 
of the system. 



3. Control parameter for avalanching systems 

We now envisage a generic avalanche model in a system of size Lq where the height 
of sand is specified on a grid, with nodes at spacing 61. Sand is added to individual 
nodes, that is, on length scale 61 at an average time rate Einj = h per node. There 
is some process, here avalanches, which then transports this dynamical quantity (the 
sand) though structures on intermediate length scales 61 < I < Lq. Sand is then lost to 
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Table 2. 11 theorem applied to an avalanching system. The sand carries a property 
with dimension S. 



Variable 


dimension 


description 




Lo 


L 


system size 




51 


L 


grid size 




e 


ST-^ 


system average 


dissipation/loss rate 


h 


ST-^ 


average driving 


rate per node 



the system (dissipated) at a time rate e over the system size Lq. The relevant variables 
for the avalanching system are given in Table 2. The two dimensionless groups are: 

ni = - = = § = f{N) (6) 

e ol 

The second parameter, 112 = f{N) is related to the number of d.o.f. of the system. The 
control parameter Hi = h/e is analogous to the Reynolds number above in that, as we 
will now show, it relates the ratio of the driving to the dissipation rates to the number 
of active, or energy containing degrees of freedom in the system. 

In Euclidean dimension D there are nodes so that conservation p2l [23t [2^ 

of the flux of sand (in an ensemble averaged sense), gives h{LQ/6l)^ ~ e which simply 
states that the rate at which sand is added to the system must on average balance 
the rate at which sand leaves. On intermediate length scales 61 < I < Lq, sand is 
transported via avalanches. There must be some detail of the internal evolution of 
the pile that maximizes the number of length scales / on which avalanches occur. For 
avalanche models this is the property that transport can only occur locally if some local 
critical gradient is exceeded; as a consequence the pile evolves through many metastable 
states. If these length scales represent d.o.f. then the number of d.o.f. available will 
be bounded by Lq and 61 so that: 

N ~ {Lo/6ir (7) 
with D > a > for > 1 {a may be fractional). We then have: 

This is in contrast to fluid turbulence since the number of d.o.f. decreases with increasing 
drive, that is, increasing = h/e. Thus we recover the SDIDT limit for SOC, namely 
Ra 0, but now explicitly identify this limit with maximizing the number of d.o.f. 
available, that is, the disorder of the system. Our result from dimensional analysis is 
to obtain Ra ~ N^'^ and to show quite generally that /Sn < 0. Following the above 
discussion of turbulence, we can go further and make the analogy Ra = Re, that is, the 
system's 'effective Reynolds number' increases with the energy/sand taken up by the 
system, i.e. with h. 

The property that the system generates many coupled d.o.f. is, for SOC, captured 
by avalanching phenomenology. This sets conditions on the microscopic details of the 
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system; specifically, there must be a separation of timescales in that the relaxation time 
for the avalanches must be short compared to the time taken for the drive to on average 
cause a cell to be come unstable so that avalanching is the dominant mode of transport. 
The critical gradient can be a random variable but provided it has a defined average value 
we have an average number of timesteps to drive a cell unstable {g5l) / {h6t) where 6t 
is the timestep. This gives two conditions for avalanching to dominate transport |25]: 




The first condition is that avalanches only occur after many grains of sand have been 
added to any given cell in the pile and is the strict SD1DT[TT1 [121 [IS] limit. However, 
if the system has large bandwidth Lq ^ 61, one can consider an intermediate behavior 
gLo ^ h6t > g6l where the driver is large enough to swamp the smallest avalanches, but 
larger avalanches persist [25]. For fixed Lq and 61, increasing h6t above g6l successively 
erodes the available d.o.f since each addition of sand swamps h6t/{g6l) cells of the pile. 
Ultimately as h and hence Ra is increased to the point where h6t ~ g6l {Lq/61)^ there 
will be a crossover to laminar flow. 

We now show that this intermediate, finite Ra behavior will be 'SOC like', with 
power law avalanche statistics sharing the same exponent as at the SDIDT limit; 
confirming our assumption above that (3^ is independent of the control parameter Ra- 
To see this, consider passing through the regime of h6t ~ g6l with h6t <C g6l {Lq/61)^ , 
which can be achieved by increasing both h and Lq such that h — > Ah and Lq 
LqA'^^/^\ This is equivalent to coarse graining the pile spatially, so that provided the 
system has self similar spatial scaling we can anticipate obtaining the same solution 
(subject to a rescaling) provided Lq LqA'^^/^\ Under nonlocal feeding and non 
overlapping avalanches [A times as many grains added at well separated positions over 
the pile) this coarse-graining may not occur. 

We illustrate this in Figures 1 and 2 with simulations of the BTW[2] sandpile in 
2D, where the driving occurs randomly in time and is spatially restricted to the 'top' of 
the pile. In all cases the critical gradient (threshold for avalanching) is (7 = 4. Figure 1 
shows two simulations of size Lq = 100 with h = [4, 16] {6t = 1 in the simulations) i.e. 
just at, and above, the regime h6t ~ g6l, but in both cases with h6t <^ g6l {Lq/61)^ . At 
/i = 16 we see that the power law statistics of the smallest avalanches is lost but there is 
still scaling over a more restricted range of avalanche sizes, i.e. we have the same scaling, 
and same exponent, but over a reduced number of d.o.f. Rescaling the avalanche sizes 
of the h = 16 run S" — > 5/16 (that is, lengthscales / I /A) recovers the behaviour of 
the h = 4 run except at the largest decade. To recover the full range we repeat the 
h = 16 run in a larger box, Lq = 400 (that is, Lq 4Lo) which is shown alongside 
the h = A, Lq = 100 run in Figure 2. Rescaling the Lq = 400 run with S S'/16, i.e. 
lengthscales 1-^1/4 then reproduces the h = 4, Lq = 100 results. This establishes 
a general property of avalanching systems that has been seen in several representative 
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SOC models[27l[25l[26l[28], see also [291 [30]. 
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Figure 1. Avalanche size normalized distributions for two runs of the 2D BTWpl [21] 
sandpile driven at the top corner formed by two adjacent closed boundaries, the other 
boundaries are open. Lq — 100 and ft, = 4 (•) and h = 16 (x); (a-left) probability 
densities; (b-right) as (a) with probability density for the h — 16 avalanche sizes 
rescaled 5 5/16. 
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Figure 2. Avalanche size normalized distributions for Lq — 100, /i = 4 (•) and 
Lq — 400, h — 16 (+); (a-left) probability densities; (b-right) as (a) with probability 
density for the h — 16 avalanche sizes rescaled S S/16. 



The above assumes the BTW case where avalanches relax instantaneously. One 
can develop this idea to introduce "running sandpiles" (as extensively studied by, for 
example, [29l ETJ EH [30]) where redistribution is no longer instantaneous. Instead, 
wherever the critical gradient is exceeded locally, / grains are moved in time 6tf. Thus 
there is a local redistribution rate per cell hf = f/Stf which we can compare with the 
driving rate h. This introduces a new dimensionless parameter [hStf)/{hf5t) which 
modifies h in the discussion above; in other words the SDIDT limit is approached 
for both /i — > and 6tf/6t — > [HI [12]. However, consistent with studies of running 
sandpiles [291 1211 [211 [30] , we can anticipate that avalanching phenomenology will persist 
for a range of finite 6tf/6t. 
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Table 3. 11 theorem applied to a simple model for an ecosystem in a space with 
Euclidean dimension D. Interactions between species processes a quantity (biomass, 
here) with physical dimension B. 



Variable 


dimension 


description 


Lo 


L 


system size 


Lc 


L 


normalization length scale 


Mp 


BT-^ 


top predator rate of consumption 






of biomass over system 


Mf 


BT-^L-^ 


rate of supply of biomass/unit volume 



Depending on the details, some SOC systems may show scaling in systems where 
the drive is in fact highly variable. One could argue (see also [TT]) that such robustness 
against fluctuations in the driver is necessary for SOC to provide a 'working model' in 
real physical systems where the idealized SDIDT limit may not be realized. 

4. Control parameter for model ecosystem 

Finally, we apply the above framework to simple models for ecosystems. We consider a 
large number of connected 'meta- species' (or groups of species/ variations occupying a 
given nichejH]) with diverse sizes and rates of predation. Each meta- species is a d.o.f. 
in the model. We insist that the details are unimportant except that each of the N 
meta- species, by acting as predator of one set of neighbors in the food web and prey to 
another set, processes some dynamical quantity, say, biomass or energy. The ecosystem 
then has a driving rate, or rate of supply H of biomass/energy per unit volume at the 
'bottom' of the web and a dissipation rate, or rate of consumption P of biomass by the 
top predators. We consider a steady state on the average which includes secular change 
in these parameters that is slow compared to the timescale for information to propagate 
through the web. For a given habitat, the abundance of species (i.e the relative number 
of distinct meta- species) grows with the size of the habitat. Although the details may 
vary, a good working approximation for the 'species- area relationship' ^32j is a power 
law, so that the number of species iV in a habitat of size Lq is given by iV ~ (-^^o)''^- ^ 
dimensionally balanced expression in D Euclidean dimensions is: 

N ~ (Lo/Lj^^ (10) 

with 7 > 0. The length scale L^. captures details of the sampling, as well as specifics of 
a given habitat and terrain. The relevant system variables are shown in Table 3. There 
are two dimensionless groups: 

= = ^^' = r = ^^^^ = ^^^^ 

thus we identify a control parameter Hi = Rb for the simple ecosystem. To relate 
this to the abundance of species we require some conservation property and to insist on 
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steady state. One possibility is to conserve some fraction of the biomass flux propagated 
through the web so that for a steady state for the system as a whole, the rate of supply 
of biomass is balanced by the rate of removal by the top predators giving L^H ~ P. 
In a system with losses, provided a fraction a of the propagated quantity is on average 
passed from one d.o.f. or species to the next, this expression is Aa^ LqH ~ P; the 
factor A also includes any recychng of the top predator biomass to the bottom of the 
web. 

We can however work with any quantity which is transferred from one species to 
another with some conservation. If we instead consider P and H to refer to integrated 
energy consumption of the top predator population and the energy taken up by the 
web per unit volume (the productivity) respectively, conservation is then the original 
'energetic- equivalence rule'|31j- that the total energy flux of a population is invariant 
with respect to body size. The control parameter Rb increases with a measure of the 
rate at which biomass (or energy) is utilized by the system as a whole (P),that is, is 
ultimately consumed by the top predator. Equivalently, it increases with the biomass (or 
energy) rate of supply to the system via the organisms at the bottom of the web, HLq\ 
these both represent the rate at which biomass/energy is processed by the ecosystem as 
a whole. We then have: 



so we obtain that Rb ~ A^^^ and quite generally that (3^ > 0. The abundance of species 
simply increases with Rb capturing the observation that diversity grows with the global 
flux of energy /biomass, that is, productivity times area [31], [32]. This result holds even 
if the species- area relationship is not a power law, it simply requires that the number 
of species grows with habitat size; intriguingly, a non- power law Pn suggests a length 
scale dependence in the abundance of species that is intermittent. 

The power law dependence implied by a power law species-area relationship suggests 
that the dependence of Rb on is rather nonlinear. A consequence is that, if we 
consider slowly increasing this control parameter in a manner that does not violate 
our assumption of a steady state, we would expect, starting from an initial state of 
few species, to see a sudden 'explosion' in diversity at some critical value Rl. This 
will depend on the details through the non universal parameter Lc] but since Lc can 
be determined through species- area abundance relationships, Rl can in principle be 
determined. 

We can consider the analogy Rb = Re, that is, we identify Rb as the ecosystem's 
'effective Reynolds number' which increases with the energy/biomass taken up by the 
system. We have then shown that, the disorder, or complexity of the ecosystem as 
expressed by the abundance of species increases with effective Reynolds number in the 
same sense as turbulence. This analogy to turbulence may be instructive in that there 
is some non universal value of the Reynolds number at which a given system makes the 
transition to turbulence. One can speculate that dynamical systems routes to turbulence 
(specifically, the Ruelle-Takens or Feigenbaum scenarios, e.g. [33]) suggest a new 
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approach to modelling the onset of the diversity of species. Understanding ecosystems 
in the context of simple models for turbulence may also provide a basis for modelling 
the "bursty" dynamics and scaling intrinsic to some ecosystems [M]. Our approach to a 
'generalized Reynolds number' outlined here potentially has wider application: to living 
organisms and societal organizations, insofar as they can be modelled |8j as webs of many 
interacting elements that process some dynamical quantity. 
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